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^ ABSTRACT 

o 

i-{ ■ We have calculated a grid of photospheric phase atmospheres of Type la super- 

^ ■ novae (SNe la) with nietallicities from ten times to one thirtieth the solar metallicity 

in the C+0 layer of the deflagration model, W7. We have modeled the spectra using 



the multi-purpose NLTE model-atmosphere and spectrum-synthesis code, PHOENIX. We 
show models for the epochs 7, 10, 15, 20, and 35 days after explosion. When compared 
^ \ to observed spectra obtained at the approximately corresponding epochs these synthetic 

spectra fit reasonably well. The spectra show variation in the overall level of the UV 
continuum with lower fluxes for models with higher metallicity in the unburned C-l-0 
layer. This is consistent with the classical surface cooling and line blocking effect due to 
metals in the outer layers of C-l-0. The UV features also move consistently to the blue 
with higher metallicity, demonstrating that they are forming at shallower and faster 
layers in the atmosphere. The potentially most useful effect is the blueward movement 
of the Si II feature at 6150 A with increasing C-l-0 layer metallicity. We also demon- 
strate the more complex effects of metallicity variations by modifying the ^^Fe content 
of the incomplete burning zone in W7 at maximum light. We briefly address some 
shortcomings of the W7 model when compared to observations. Finally, we identify 
that the split in the Ca H-l-K feature produced in W7 and observed in some SNe la is 
due to a blending effect of Ca II and Si II and does not necessarily represent a complex 
abundance or ionization effect in Ca II. 



Subject headings: line: formation — nuclear reactions, nucleosynthesis, abundances 
radiative transfer — stars: atmospheres — supernovae: general 



1. Introduction 

Type la supernovae (SNe la) are recognized for their near uniformity as standard candles. This 
has led to the use of SNe la in cosmology (Branch 1999; Garnavich et al. 1998; Riess et al. 1998; 
Perlmutter et al. 1999). The development of empirical calibrations between peak brightness and 
light curve shape (Phillips 1993; Hamuy et al. 1996; Riess et al. 1996; Perlmutter et al. 1997a,b; 
Tripp (fe Branch 1999), have enhanced the usefulness of SNe la as distance indicators. Nevertheless 
theoretical modelers have yet to agree on the source of these apparently systematic variations. A 
primary concern is the evolutionary lifetime of SNe la progenitors and the possibility of significant 
deviations of distant SNe la from their well observed counterparts in the local galactic neighborhood. 
If galactic chemical evolution occurs slowly, then more distant SNe arise from a younger, metal poor 
population. On the other hand, the metallicity variations in the local sample may already span the 
range of the entire observational sample. 

We probe the possible effects of progenitor metallicity variations on the observed spectra, 
by modifying the parameterized deflagration model, W7 (Nomoto et al. 1984; Thielemann et al. 
1986). Using base fits to observations (§ 2.3) we have scaled all elements heavier than oxygen in 
the unburned C+0 layer of W7 to simulate the effects of various metallicities in the progenitor 
system. Hoflich et al. (1998) explored this question by modifying the pre-explosion metal content of 
a particular SNe la model and noted the differences in final composition, light curves, and spectra. 
They found that the composition of the partially burned layers of the ejecta yielded larger quantities 
of ^^Fe. A similar effect can be seen in the lowered ^^Fe abundance when W7 is calculated using 
a pure C+0 mixture without other metals (Iwamoto et al. 1999). We have modified the ^^Fe 
abundance of the partially burned layers of W7, to replicate this effect. While Hoflich et al. (1998) 
focused mainly on the effects of metallicity variations on the light curve and energetics, here we 
concentrate exclusively on its effects on the observed spectra, particularly at early times where the 
formation of the spectrum occurs in the unburned C+0 layer, which is most sensitive to initial 
progenitor metallicity. 

Our computational methods are given in § 2. In § 3 we show the effects of metallicity modifi- 
cation of the C+0 layer on the synthetic spectra at each epoch. In § 4 we show the effects of ^^Fe 
on day 20 and 35 model spectra. In § 5.1 we discuss the evolution of the Si II feature at 6150 A. 



2. Methods 

2.1. Abundance Modification 

To model the effects of metalhcity on SNe la, we have modified the base W7 model at several 
epochs. W7, while clearly not the complete model for SNe la, is a good starting point for spectral 
modeling. The composition structure of W7 as a function of velocity does reasonably well in repro- 
ducing the observed spectra of SNe la (Harkness 1991a, b; Branch et al. 1985; Nugent et al. 1997). 
By making separate calculations to consider the direct (C+O metalhcity) and the indirect explosive 
nucleosynthetic consequences (^^Fe enhancement), our computational methods allow us to probe 
these effects of progenitor metalhcity independently. Hoflich et al. (1998) changed the composition 
before computing the hydrodynamics and explosive nucleosynthesis of the model. This gives them 
consistent nucleosynthesis and energetics within the context of their chosen input physics. How- 
ever, it is difficult to separate effects in spectra that arise from different consequences of the initial 
metalhcity variation. 

To modify the metalhcity of the unburned C-l-0 layer we have scaled the number abundances 
of elements heavier than oxygen in the velocity range 14800-30000 km s~^ by a constant factor (^, 
such that for all species i heavier than oxygen, the new number abundance n[, is given by 

n'i = Cui. 
The mass fractions are then renormalized, 

in each layer. For all epochs we have used C-l-0 metalhcity factors, (", of 10, 3, 1, 1/3, 1/10, and 
1/30. These effects are discussed in detail in § 3. 

The nucleosynthesis results of Hoflich et al. (1998) show that modifying the metalhcity of the 
progenitor white dwarf changes the quantity of ^^Fe produced in the incomplete burning zone in 
a manner approximately proportional to the metalhcity. This is due to excess neutrons in the 
progenitor over a pure C-l-0 mixture, particularly ^^Ne (see for example Arnett 1996; Nomoto 
et al. 1984; Iwamoto et al. 1999, and references therein). To simulate this effect in W7, we have 
scaled the ^^Fe abundance in the incomplete burning zone of W7, 8800-14800 km s^^, in the same 
manner as the C-|-0 layer metalhcity. 

The mass fractions are then renormalized in each layer. We have used ^^Fe abundance factors, ^, of 
3, 1, 1/3, 1/10, and 1/30 for the epochs 20 and 35 days after explosion. For the remaining epochs 
we have only computed models for 1/10 ^^Fe abundance (^ = 0.1). These effects are discussed in 
detail in § 4. 

Finally, we simulate the full effects of progenitor metalhcity on the output spectra, by com- 
bining the two effects, C+O layer metalhcity and incomplete burning zone ^^Fe abundance, into 



a single series of models. For each model we have used the same factor for both C+0 metallicity 
and ^^Fe abundance (^ = (). These factors are the same as for the ^'^Fe abundance modifications 
above. These effects are discussed in detail in S 4. 



2.2. Computational Methods of PHOENIX 

The calculations were performed using the multi-purpose stellar atmospheres program PHOENIX 9 . 1 
(Hauschildt & Baron 1999; Baron & Hauschildt 1998; Hauschildt et al. 1997a,b, 1996). PHOENIX 
solves the radiative transfer equation along characteristic rays in spherical symmetry including all 
special relativistic effects. The non-LTE (NLTE) rate equations for several ionization states are 
solved including the effects of ionization from non-thermal electrons from the 7-decay energy of 
the ^^Ni core. The atoms and ions calculated in NLTE are: He I (11 levels), He II (10), C I (228), 
O I (36), Na I (3), Mg H (18), Si H (94), S H (85), Ti H (204), Fe H (617), and Co H (255). Each 
model atom includes primary NLTE transitions, which are used to calculate the level populations 
and opacity, and weaker secondary NLTE transitions which are are included in the opacity and 
implicitly affect the rate equations via their effect on the solution to the transport equation. In 
addition to the NLTE transitions, a number of LTE line opacities for atomic species not treated 
in NLTE are treated with the equivalent two-level atom source function, using a thermalization 
parameter, a = 0.1, as in Nugent et al. (1997) for LTE lines for all models in this paper. The 
atmospheres are iterated to energy balance in the co-moving frame, 
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where, S is the deposited energy from radioactive decay, 7 and f3 have their usual meaning and 
c = 1. This equation neglects the explicit effects of time dependence in the radiation transport 
equation, however the term on the right hand side implicitly includes these effects. The models are 
parameterized only by the day, which determine the radii and amount of radioactive decay, and by 
the luminosity parameter, rj, which is defined as, 
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where L*^^ is the instantaneous deposition of radioactive decay 7-ray energy. The deposition 
of 7-rays is determined by solving a grey transport equation (Nugent 1997; Nugent et al. 1997; 
Sutherland & Wheeler 1984) as a function of time. This approach has been shown to be of sufficient 



accuracy for our purposes when compared to detailed Monte-Carlo calculations (Swartz et al. 1995). 
We use two different inner boundary conditions for our models. For days 7, 10, and 15, which are 
optically thick at the core, we use a diffusive inner boundary condition. For days 20 and 35 which 
have lower total optical depth we use a 'nebular' inner boundary. This boundary condition takes 
the downward intensity of any ray and assumes that it passes unimpeded through the core of the 
atmosphere, becoming the upward intensity on the other side, with only the necessary Lorentz 
transformation. Additional details on the use of PHOENIX in modeling SNe la atmospheres can be 
found in Nugent (1997). 



2.3. Baseline Models 

To perform numerical experiments on the metallicity effects in the spectra of SNe la, we 
need reasonable base models for each epoch. For each day, an observed spectrum of SN 1994D 
was fitted with a W7 model of appropriate age. Several models were computed with various 
luminosity parameters, rj, to fine tune the spectral shape. A separate paper fitting W7 to SN 1994D 
photospheric era spectra (Lentz et al. in preparation) shows the comparisons to observations. 
Hoflich &: Khokhlov (1996) find rj (their Q) in the range 0.7 < rj < 1.8. When rj > 1 stored 
radiative energy is being released and when r/ < 1 radiative energy is being stored. 

The maximum light (epoch day for the observations) fit for SN 1994D to W7 (with rj = 1.0) 
at day 20 after the explosion is the same as in Nugent et al. (1997). For the two model fits at 
days 10 (which we compare to day -9 with respect to B maximum) and 15 (compared to day -4 
with respect to B maximum) for W7 we have used luminosity parameters of r/ = 0.4 and rj = 0.8 
respectively. For these two models the fits are generally good, but with all of the pre-maximum 
spectra, the red edge of the Si II feature at 6150 A does not extend far enough to the red. For 
model day 7 (fit to day -12) we chose the the model with luminosity rj = 0.1. This model fits the 
luminosity and most features well; however, like the day 10 models the Ca H-l-K feature not extend 
blueward enough when compared to the earliest observations. Hatano et al. (1999) have modeled 
the same spectra and found that the Ca H-l-K feature requires ionized calcium with velocities up 
to ~ 40000 km s~^. We have conducted experiments extending the density profile of W7. This 
did not provide enough optical depth in the Ca H-|-K feature to affect the line profile. Neither do 
our models reproduce the secondary feature at 4700 A that Hatano et al. (1999) ascribe to Fe II 
absorption in the C-|-0 layer. For our post-maximum model, day 35 (15 days after maximum), we 
fit the observation with a model with luminosity, r/ = 1.5. This model fits generally; however, the 
large feature at 4800 A (probably Fe II) is not seen in the observations. The observed Na I D-line 
feature is missing, probably due to deficiencies in the sodium model atom, which has already been 
improved in the latest version of PHOENIX. Here, we are interested in small differential effects so it 
was essential that all models be run with the code version frozen. 



3. Metallicity of Unburned C+O Layer 

The simplest effect of metallicity on the spectral formation in SNe la is the change in the metal 
content of the unburned C+O layer neglecting changes in the density structure and deeper layers. 
In this section we examine these effects independently of other effects of progenitor metallicity on 
SNe la. Since ongoing supernova searches are expected to discover SNe early, and early spectra probe 
the outermost layers only, this approach is sensible and yields physical insight that is somewhat 
model independent. 

When we look at the overall UVOIR synthetic spectra (Figures 1, 6, 8, 11, and 13) of the 
models with variations in the C+O layer metallicity, we see two consistent and significant effects: 
shifts in the UV pseudo-continuum level (expanded view for day 7 in Figure 2) and variations in 
the Si II line at 6150 A (expanded views in Figure 5). 

The general effect in the UV is the increase in the UV pseudo-continuum level with decreasing 
metallicity. Simultaneous is the redward (blueward) shift of most UV features with decreasing 
(increasing) metallicity. In the UV, the line-forming region is in the C+O layer. As the metallicity 
decreases, the line forming region must reach deeper into the atmosphere to have the same line 
opacity, resulting in smaller line velocities. Modification of the C+O layer metal abundance gives 
a classic surface cooling effect, lower temperatures for higher metallicity. The higher temperatures 
of the lower metallicity C+O atmospheres give higher thermal fluxes, moving the UV pseudo- 
continuum higher with lower metallicity. The surface cooling and the resulting shifts in UV pseudo- 
continuum are evident at every epoch. There is the complementary effect of additional metals 
increasing the line blocking. We make no attempt to separate these two effects in this paper. 

The Si II line at 6150 A (Figure 5) shifts blueward with increases in metallicity for epochs 
through day 20. These shifts demonstrate that some line formation in this feature takes place in 
the C+O layer. The earlier epochs show large variations in the total depth of the feature which 
implies that the line forms less in the incomplete burning zone with its large, unchanging silicon 
abundance, and more in the C+O layer where the silicon abundance changes. At later epochs these 
conditions are reversed, resulting in smaller changes with C+O layer metallicity variation. These 
effects are discussed further in § 5.1. 

The Mg II "h+k" feature at 2600 A (Figure 2) does not move to the blue or red as the metallicity 
varies. The decrease in Mg II h+k feature strength with increasing metallicity is caused by the 
increasing UV line blanketing from background line opacity in the C+O layer. The Mg II absorption 
occurs mostly in the deeper, partially burned layer that is highly enriched in magnesium. We have 
confirmed this hypothesis by calculating diagnostic output spectra without using any background 
opacity (see Baron et al. 1999, for a discussion of the method). These diagnostic spectra suggest 
that this feature (and the feature near 2800 A) may actually be due to a complicated blend of 
Mg II UV 1-4 transitions: h+k (UV 1) A2798, (UV 2) A2796, (UV 3) A2933, and (UV 4) A2660. 



3.1. Day 7 

Our grid of synthetic spectra of W7 on day 7 with C+0 metahicity variations are shown 
in Figure 1. The case (" = 10 is hkely extreme and is shown for illustrative purposes only, we 
don't consider it further at this epoch. The UV spectra (Figure 2), show variations in the pseudo- 
continuum due to surface cooling, line blocking, and feature shifts due to changes in the feature 
formation depth (velocity) . The pseudo-continuum level in the optical and IR vary with metallicity. 
Lower metallicities have lower pseudo-continuum fluxes. Figures 3 and 4 display optical and infrared 
spectra which illustrate this effect more clearly. This is due to backwarming. Figure 3 shows that 
the line features at 3650 A (Ca II H-^K), 4250 A (Fe II, Mg II, k Si III), 4700 A (Fe II), and 
5650 A (Si II) nearly disappear as the metallicity drops to 1/30 of normal. The 4700 A Fe II 
feature is weak and doesn't extend far enough to the blue in the base model fits (§ 2.3). Increasing 
the metallicity by even a factor of 10 does not extend the feature blueward enough (to sufficiently 
large velocities) to match the observations. The same is true of the Ca II H-l-K lines. This indicates 
that more mass extending to higher velocities is needed to fit this feature. The Si II feature at 
6150 A (Figure 5a) shows the blueward movement of the feature minimum and blue-edge wall 
with increasing metallicity. The slope of the red edge of the feature changes and at the lowest 
metallicities the feature is not strong enough to form the red emission wing of the P Cygni feature. 
The infrared spectra (Figure 4) show several features that weaken with lower metallicity, Mg II 
A9226 A and Si I complexes at A10482 and A10869 A (cf. Millard et al. 1999), as well as two features 
that strengthen with lower metallicity. These two O I features (7400 A and 8150 A) weaken as the 
oxygen abundances decreases in favor of higher Z metals. 



3.2. Day 10 

Our grid of synthetic spectra of W7 on day 10 with metallicity variations in the C-l-0 layer 
(Figure 6) display similar variations to those in the day 7 spectra (Figure 1). The UV variations 
again show the blue-shifted features and lower continua with increasing metallicity. The optical 
pseudo-continuum effects (Figure 7) are still present, but smaller. We again see lines that deepen 
with increased metallicity such as 3650 A (Ca II H-l-K), 4250 A (Fe II, Mg II, & Si III), and 6100 A 
(Si II), but these effects are less dramatic than at day 7. The Si II feature (Figure 5b) still shows 
the changing depth and slopes of the line edges with metallicity, but, at day 10 the depth of the 
Si II feature in the lowest metallicity model is larger than in the same model on day 7, relative 
to the highest C-l-0 layer metallicity in the respective epoch. These effects on the optical lines 
indicate that line formation for certain strong lines is now taking place in layers below the C-l-0 
layer. The effects of line formation below the C-l-0 layer becomes more important as the supernova 
atmosphere expands and becomes less opaque. 



3.3. Day 15 

Our grid of synthetic spectra of W7 on day 15 with C+0 metahicity variations (Figure 8) 
are again similar to those in the day 7 and 10 spectra (Figures 1 & 6). The decreasing effects of 
C+0 layer metallicity are apparent. The UV pseudo-continuum variation with metallicity remains 
strong. The optical (Figure 9) and the near infrared (Figure 10) show the backwarming optical/IR 
pseudo-continuum flux effect which occurred in the optical at day 10, and extended well into the IR 
at day 7. Features which increased in strength with increasing C+0 metallicity are a blend of Fe II, 
Mg II, & Si III at 4250 A, the multi-species blend at 3300 A, Mg II at 8850 A, and an unidentified 
feature at 10400 A (possibly Si I). The feature at 3300 A is a blend of weak lines that forms in the 
C+0 layer. Figure 9 illustrates the blue shifting of features as changes in metal content move the 
depth (and thus velocity) of line formation. 

The two component feature at 3700-3900 A is usually labeled Ca II H+K. Nugent et al. (1997) 
have identified this 'Ca split' as arising from a blend of Ca II H+K with Si II A3858. We have 
calculated a series of diagnostic spectra at this epoch using the same temperature structure, but 
without background opacity. We have confirmed the identification of the blue wing with the Si II 
line. When spectra are calculated without any other line opacities, Ca II forms a pair of features. 
The blue Ca II absorption forms from the A3727 line, the red absorption from the H+K lines. The 
Si II feature falls on the peak between the two Ca II absorptions. As the Si II feature strengthens, 
the flux displaced by its absorption 'fills-in' the Ca II H+K absorption creating a 'split'. With 
enough Si II only the blue absorption feature remains, while the red feature becomes an inflection. 
The lack of a split in some SNe la indicates that either the Si II feature is weaker or the Ca II 
H+K is stronger, preventing the formation of a split. This helps to illustrate that in supernovae of 
all types, features are often the result of more than one multiplet, or even ionic species. 

The Si II feature at 6150 A (Figure 5c) shows much smaller, but still significant, effects due to 
C+0 metallicity. The blue edges are parallel and the red edges nearly so. The depth contrast is 
now much smaller than before. This is strong evidence that the feature is forming primarily below 
the C+0 layer at this epoch. Figure 10 shows the displacement of oxygen by metals in the O I 
lines at 7450 A, 8100 A, and 8250 A, the latter two of which are superimposed on the stronger 
Ca II 'IR-triplet' which, like the H+K absorption, forms in the calcium rich incomplete burning 
zone and is unaffected by the C+0 metallicity. 



3.4. Day 20 

Our grid of synthetic spectra of W7 on day 20, approximately maximum B magnitude, with 
C+0 metallicity variations (Figure 11) shows the continuing reduction in importance of the C+0 
layer metallicity as the atmosphere becomes more optically thin. The UV pseudo-continuum and 
feature shift effects are still strong. The only remaining features which display C+0 layer metallicity 
effects in the optical (Figure 12) are the Ca II H+K, the 4250 A (Fe II, Mg II, & Si HI), and 



6150 A Si II features. As at day 15, the 'Ca split' remains, the blue component increasing in 
strength and the red component decreasing with increasing C+0 metallicity. The Si II feature 
at 6150 A (Figure 5d) now shows parallel edges (both blue and red) for the various C+0 layer 
metallicities. The relatively small changes in the line strength indicate that the feature forms 
mostly in the deeper, silicon rich layers, but some measurable effects due to C+0 layer metallicity 
are still evident. 



3.5. Day 35 

Our grid of synthetic spectra of W7 on day 35 with C+0 metallicity variations are plotted in 
Figure 13. The UV features still show the same pseudo-continuum and line shifting behavior seen 
in the previous epochs. This illustrates that even at this epoch the pseudo-continuum formation in 
the UV is still in the unburned C+0 layer. The Si II feature at 6150 A has only small variations 
which can not be separated from the pseudo-continuum effects. It should be noted that the Si II 
feature at this epoch fits the observations rather poorly. 

In the optical we see the blue-shifting of the 4800 A (Fe II) feature with increasing metallicity 
without significant changes to the overall width, shape, or depth of the feature. This indicates 
that it forms mostly at deeper layers, with smaller effects from the C+0 layer. When looking for 
a base fit for this epoch, we find that while the line shapes near 5000 A are better with somewhat 
hotter models, those models had poor overall spectral shape or color. To diagnose this sudden 
change in model behavior, we have plotted the temperature profile of the three models with the 
lowest metallicity, the density profile, and the carbon abundance (Figure 14). We can see that the 
^ = 1/30 model has a definite temperature inversion. This inversion corresponds to the position of 
the density spike in W7 arising from the deflagration wave. The higher density coincides with the 
change from the more efficient cooling available by intermediate mass elements to that of the less 
efficient cooling of C+0 creating the temperature inversion. In models with higher metal content 
the effectiveness of cooling by metals in the C+0 layer provides the needed cooling to prevent the 
formation of an inversion. 



4. ^^Fe Content of Incomplete Burning Zone 

4.1. Pre-Maximum Light Epochs 

We have computed models with ^^Fe abundance reductions of 1/10 in the intermediate burning 
zone and models with the combined 1/10 C+0 metallicity and 1/10 ^^Fe abundance reductions 
for the epochs 7, 10, and 15 days. When we compared the ^^Fe reduced abundance models with 
the related models containing the same C+0 metallicities (normal and 1/10 normal, respectively) 
we found no changes in feature strength. Some slight differences in pseudo-continuum levels were 
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seen in the day 10 and 15 comparisons. These may be due to weaker versions of the ^^Fe abun- 
dance caused backwarming described in § 4.2. Since each comparison pair includes two models 
with different composition in the incomplete burning zone this may slightly affect the temperature 
structure. 



4.2. Day 20 

Our grid of synthetic spectra of W7 on day 20, approximately maximum B magnitude, with 
incomplete burning zone ^^Fe abundance variations is shown in Figure 15. Some small vertical 
displacements of the UV flux, without blue- or red-shift, can be understood by the effects of surface 
cooling/backwarming in the incomplete burning zone. As intermediate mass elements are replaced 
by ^^Fe, the larger line opacity of iron cools the ^^Fe rich models. Since this does not affect the 
temperature gradient in the C+0 layer, the temperature shift remains constant throughout the 
C+0 layer. 

In the optical (Figure 16) we can see a few features that decrease with increasing Fe, such as 
the Si II and S II in the 4850 A feature, the S II "W" at 5300 A, and the Si II feature at 6150 A. 
This is caused by displacement of the species forming the line by additional ^^Fe. Several features 
can be seen to strengthen with greater ^^Fe abundance. This indicates that they are formed at 
least in part by iron in the incomplete burning zone. These features included the feature at 4100 A 
that erodes the peak of a neighboring feature and the red wing of the 4300 A feature. 

Our grid of synthetic spectra of W7 on day 20 with incomplete burning zone ^^Fe abundance 
variations and C+0 layer metallicity are shown in Figure 17. The effects of the combined modi- 
fications mostly separate into the effects seen in Figures 11 &: 15. The UV displacements in flux 
and the effects on the Si II 6150 A feature in the combined models reflect the effects of metallicity 
variation in the C+0 layer alone. The S II "W" at 5300 A and the "peak erosion" line at 4100 A 
show the primary effect of ^^Fe abundance on the optical spectra (Figure 18) when both the ^^Fe 
abundance and the C+0 layer metallicity are varied simultaneously. The effect of the combined 
modifications on the Ca II H+K feature are small. The 4350 A feature in the combined modifica- 
tion has the combination effects of deeper red wing strength with increasing ^^Fe abundance and 
a deeper blue wing (Mg II) with increasing C+0 layer metallicity. The changes in the C+0 layer 
metallicity and in the ^^Fe abundance of the incomplete burning zone each have effects which are 
separate from one another and combined effect is essentially the sum of the two effects. 



4.3. Day 35 

Our grid of synthetic spectra of W7 on day 35 with incomplete burning zone Fe abundance 
variations are shown in Figure 19. The UV pseudo-continuum varies due to the 'surface' cooling 
and additional line blocking in the incomplete burning zone. The O I feature at 8000 A and the 
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Si II feature at 6150 A become weaker as oxygen and silicon are displaced by ^^Fe. There are 
several significant effects in the infrared, however, observations which extend far enough into the 
IR to compare with real SNe la are not available. 

Our grid of synthetic spectra of W7 on day 35 with combined incomplete burning zone ^ Fe 
abundance and C+0 layer metallicity variations are shown in Figure 20. The effects of the combined 
modifications mostly separate in the the effects shown in Figures 13 &: 19. The UV pseudo- 
continuum variation, changes near the 5000 A Fe II feature, and the changes to the llOOOA feature 
reflect the C+O metallicity modification alone. The small shift in the Si II feature at 6150 A and 
the surrounding pseudo-continuum show the contributions of ^^Fe abundance modification. 



5. Discussion 

5.1. Evolution of Si II 

Figure 5 shows the evolution of the Si II 6150 A feature to maximum light. The feature grows 
stronger and steeper as line formation of the Si II A6355 feature moves into the silicon rich layers 
of W7. As line formation moves into the silicon rich zones, the effects of C-l-0 layer metallicity on 
line formation, are reduced, but not eliminated. The blue-shift velocities of the deepest points of 
the Si II feature are plotted for these models in Figure 21. Except for the extreme case, (" = 10, 
the blue-shift velocities increase monotonically with C+O layer metallicity through day 20. The 
increasing opacity from the C+O layer moves the feature blueward. The velocity shifts due to 
metallicity are degenerate with changes that could be expected from silicon rich material extending 
to higher velocities. However, the effects of primary line formation in the C+O versus silicon rich 
layers can be distinguished by line shapes. These general trends and spreads in blueshift velocities 
are similar to those seen in the data by Branch & van den Bergh (1993). They found that the 
slower blueshift velocities tended to be found in earlier galaxy types. A similar study correlating 
blueshift velocities, peak magnitude (or a suitable proxy), and a more quantitative estimate of the 
pre-supernova environment metallicity will be the subject of future work. 



5.2. Other Features and Eflfects 

The UV pseudo-continuum shows the effects of metallicity through the surface cooling of the 
C+O layer, additional line blocking and the shifting of lines which form at faster moving layers with 
higher metallicity. The UV pseudo-continuum displacement is relatively constant and still present 
at 35 days after explosion. The UV displacement over the entire range of models is typically ~ 0.5 
dex. For the more likely range of metallicities, 1/3 to 3 times solar, the change is up to ~ 0.2 
dex or 0.5 magnitude. Unfortunately, we cannot use near-UV flux as a metallicity indicator. The 
UV flux is diagnostic of the temperature, density, and radius (velocity) of the C+O layer, but is 
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not uniquely determined for any one quantity. The related backwarniing causes pseudo-continuum 
shifts in the optical and near infra-red for the early epoch spectra. This backwarming shift in 
continua fades in strength as pseudo-continuum formation moves to deeper layers. A smaller surface 
cooling/backwarming effect exists in the partially burned layers from changes in ^^Fe abundance, 
beginning at day 20. 

Hoflich et al. (1998) also report a change in UV pseudo-continuum for models with different 
metallicity by noting decreases in the U band magnitude with increasing metallicity. They show a 
change in flux in the UV spectra presented, but the change occurs in the opposite direction with 
metalicity to that which we find here. This is due to the difference in the density structures between 
delayed detonation models (DDs) they employ and that of the parameterized deflagration model 
W7 that we use here. In the DD model the lower metallicity model forms the UV pseudo-continuum 
at a smaller radius and so the flux is lower for smaller metallicities, since the radial effect wins out 
over the opacity. 

Some features which had significant impact from C+0 metallicity at early epochs are nearly 
unaffected at later epochs. At day 7 we find that most of the optical features nearly disappear when 
the metal content is dropped to 1/30 normal. This indicates that much of line formation for these 
features takes place in the C+0 layer. In later epochs the impact of the C+0 layer metallicity on 
most optical features becomes quite small, since line formation occurs mostly in the deeper layers. 

The influence of ^'^Fe abundance in the incomplete burning zone is very small at early epochs, 
since significant spectrum formation occurs in the C+0 layer. At maximum light and later, ^^Fe 
abundance variations change Fe feature strengths, have some small temperature based effects on 
the pseudo-continuum, and change the strength of certain features from species, e.g. sulfur, which 
are displaced by ^^Fe abundance changes. 

We have confirmed the Nugent et al. (1997) identification of Si II A3858 as a component of the 
Ca II 'split'. To the blue of the Ca II H+K feature is the emission peak from another Ca II line. 
When the Si II line forms, the Si II absorption falls on that peak and the emission from the Si II 
line forms a peak in the center of the 'split'. For models with a strong Si II feature, the absorption 
minimum of the red wing of the Ca II H+K will seem to disappear. For W7, changes in metallicity 
alter the strength of the Si II line and thus, the shape of the Ca 'split' feature. In real SNe la, 
other factors may prevent the formation of a 'split' such as stronger Ca II H+K, weaker Si II, or 
temperature effects. Jha et al. (1999) show spectra for several normal SNe la. The Ca II H+K 
feature in the observations display a range of morphologies similar to those seen in our synthetic 
spectra. We do not require any abundance or ionization effect to produce this feature. While 
this interpretation is strictly correct in the model W7, we suspect it also produces the observed 
feature in SNe la, however, detailed comparisons of the calculated velocities to those of the observed 
features are required. 
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6. Conclusions 

By calculating a series of model atmospheres with abundance variations around the base W7 
model for SNe la, we have demonstrated unexpected and complex effects on the output spectra. 
The UV spectra show lower flux and blueshifting lines as surface cooling, additional line blocking, 
and outward movement of the line-forming region occur with higher metallicity in the unburned 
C+0 layer. We have demonstrated, at epochs well before maximum light that line formation 
occurs largely in the C+0 layer for species that will later form in newly synthesized material. The 
'splitting' of the Ca II H+K feature we can now better understand as a blend of a Si II feature 
with the stronger Ca II lines, without abundance or ionization effects in the W7 model. We have 
shown that the strength, profile, and velocity of the SNe la characteristic Si II 6150 A feature are 
affected by C+0 layer metallicity. This provides a mechanism for the variation in the blueshift of 
the Si II feature without variations in the explosion energy. For exploding white dwarfs of different 
metallicities Hoflich et al. (1998) and Iwamoto et al. (1999) find changes in nucleosynthesis. We 
tested these effects on spectra by varying the ^^Fe abundance in the incomplete burning zone. We 
have found that the ^'^Fe abundance has negligible effect on pre-maximum spectra and relatively 
little effect afterwards. The effects of progenitor metallicity variations can mostly be separated into 
effects due to ^^Fe, and those due to C+0 metallicity. 

Previous spectral studies (Branch et al. 1985; Harkness 1991a, b; Nugent et al. 1995, 1997) have 
shown that W7 has the appropriate composition structure to reproduce photospheric era spectra. 
In preparing the base fits for this paper we have found that the quality of the fits decreases away 
from maximum light. The temperature structure inversion in the lowest metallicity model at day 
35 demonstrates the importance of changes (w 20%) in the temperature structure. We see that 
density/temperature structures are important in fitting spectral features. In forthcoming work we 
will present the results of numerical experiments to alter the density structure of W7 to make the 
model better correspond to the observations. 

We have shown through parameterized abundance modifications of the SNe la model W7 
that pre-explosion metallicity can have detectable effects on the output spectra at every epoch. 
Due to the uncertainty in hydrodynamic models and severe blending of lines in SNe la spectra, 
we cannot give a prescriptive analysis tool for measuring the pre-explosion metallicity of SNe la. 
However, hopefully this can contribute to the understanding of the diversity of SNe la, and the ways 
that various progenitor metallicity effects can affect SNe la. Studies computing detailed spectra 
of hydrodynamic models that include progenitor metallicity in the evolution of the star and the 
hydrodynamics and nucleosynthesis of the supernova should show what effects overall progenitor 
metallicity have in creating SNe la diversity. 

This work was supported in part by NSF grants AST-9731450 and AST-9417102; NASA grant 
NAG5-3505, an IBM SUR grant to the University of Oklahoma; and by NSF grant AST-9720704, 
NASA ATP grant NAG 5-3018 and LTSA grant NAG 5-3619 to the University of Georgia. Some 
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Fig. 1. — Models with various metallicities in C+0 layer at 7 days after explosion. Thick solid 
denotes 10 times the normal C+0 metallicity, thick dotted — 3 times, short dashed — normal, long 
dashed — 1/3, dot-dashed — 1/10, and thin solid — 1/30. 
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Fig. 2. — Expansion of ultraviolet region of Figure 1 
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Fig. 3. — Expansion of optical region of Figure 1 
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Fig. 4. — Expansion of infrared region of Figure 1 
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Fig. 5. — Expansion of Si II 6150 A region of a. Figure 1 (Day 7), b. Figure 6 (Day 10), c. Figure 
(Day 15), and d. Figure 11 (Day 20). 
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Fig. 6. — Models with various metallicities in C+0 layer at 10 days after explosion. Labels are the 
same as Figure 1. 
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Fig. 7. — Expansion of optical region of Figure 6 
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Fig. 8. — Models with various metallicities in C+0 layer at 15 days after explosion. Labels are the 
same as Figure 1. 
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Fig. 9. — Expansion of blue regions of Figure 
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Fig. 10. — Expansion of near infrared region of Figure 8 



26 



14.5 



■y 13.5 



12.5 




1500 2500 3500 4500 5500 6500 7500 8500 

X (Angstroms) 



Fig. 11. — Models with various metallicities in C+0 layer at 20 days after explosion. Labels are 
the same as Figure 1. 
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Fig. 12. — Expansion of optical region of Figure 11 
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Fig. 13. — Models with various metallicities in C+0 layer at 35 days after explosion. Labels are 
the same as Figure 1. 
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Fig. 14. — Panels from top to bottom display carbon abundance, density, and temperature profiles 
for low metallicity models at day 35. The line styles for the temperature profiles are the same as 
for the corresponding spectra in Figure 13. 
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Fig. 15. — Models with various ^^Fe abundances in the incomplete burning layer at 20 days after 
explosion. Solid denotes 3 times the normal ^^Fe abundance in the incomplete burning zone, thick 
dotted — normal, short dashed 1/3, long dashed — 1/10, and dot-dashed — 1/30. 
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Fig. 16. — Expansion of optical region of Figure 15 
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Fig. 17. — Models with various ^^Fe abundances in the incomplete burning layer and metallicities 
in the C+0 layer at 20 days after explosion. These models combine the effects in the models in 
Figures 11 & 15. Solid denotes 3 times the normal ^^Fe abundance in the incomplete burning zone 

and C+0 layer metallicity, thick dotted — normal, short dashed 1/3, long dashed — 1/10, and 

dot-dashed — 1/30. 
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Fig. 18. — Expansion of optical region of Figure 17 
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Fig. 19. — Models with various ^^Fe abundances in the incomplete burning layer at 35 days after 
explosion. Labels are the same as in Figure 15. 
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Fig. 20. — Models with various Fe abundances in the incomplete burning layer and metallicities 
in the C+0 layer at 35 days after explosion. These models combine the effects in the models in 
Figures 13 & 19. Labels are the same as in Figure 17. 
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Fig. 21. — Blueshift velocities of Si II 6150 A feature as a function of epoch. Down triangle symbols 
represent 1/30 normal C+0 layer metallicity, left triangles 1/10 normal, up triangles 1/3 normal, 
diamonds normal, squares 3 times normal, and circles 10 times normal metallicity. Data are from 
lines shown in Figure 5. 



